Optimization strategies of chromosome search- and- capture by microtubules: exact 

analytical results for a single fixed target 
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The mitotic spindle is an important intermediate structure in eucaryotic cell division, in which each 
of a pair of duplicated chromosomes is attached through microtubules to centrosomal bodies located 
close to the two poles of the dividing cell. It is widely believed that the spindle starts forming by the 
'capture' of chromosome pairs, held together by kinetochores, by randomly searching microtubules. 
We present a complete analytical formulation of this problem, in the case of a single fixed target 
and for arbitrary cell size. We derive a set of Green's functions for the microtubule dynamics and 
an associated set of first passage quantities. An implicit analytical expression for the probability 
distribution of the search time is then obtained, with appropriate boundary conditions at the outer 
cell membrane. We extract the conditions of optimized search from our formalism. Our results are 
in qualitative and semi-quantitative agreement with known experimental results for different cell 
types. 

PACS numbers: 



I. INTRODUCTION 

Microtubules are one class of polymeric filaments in the 
eucaryotic cell, whose sub-unit is a hetero-dimer of alpha- 
and beta-tubulin. Microtubules therefore possess struc- 
tural polarity, and the ends are differentiated as plus and 
minus ends. A hall-mark of microtubules is their unique 
mechanism of assembly and disassembly: a polymerizing 
microtubule can abruptly start shrinking by losing sub- 
units and vice-versa, a process referred to as dynamic 
instability (reviewed in [3]). The stochastic switching pro- 
cess between growth and shrinkage is referred to as catas- 
trophe and the reverse process is called rescue. Between 
a rescue and a catastrophe, a microtubule grows in length 
by polymerizing and between a catastrophe and rescue, 
it shrinks. In vivo, a third state called pause is also 
observed where the length remains static. Microtubules 
usually nucleate from organizing centers called ccntro- 
somes, but may also be found free in the cytoplasm. 

Microtubules play a central role in eukaryotic cell di- 
vision. An important milestone in the cell division cycle 
is the formation of the metaphase spindle, where all the 
duplicated chromosome pairs, held together by kineto- 
chores are aligned along the cell 'equator (the "metaphase 
plate" ) in such a way that each chromosome of a pair is 
facing one of the poles of the cell, and attached to one or 
more microtubules emanating from a centrosome located 
near that pole. How is the spindle assembled? It is gen- 
erally believed that microtubules nucleating from each 
centrosome randomly searches the surrounding space for 
chromosomes by alternately growing and shrinking (the 
random search-and-capture model, and are stabilized 
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upon contact with a kinctochoreQ, which has also been 
supported by experimental observations Q. However, in- 
vestigations over the last decade or so have revealed that 
the chromosomes do not remain passive in this process; 
rather, the kinetochores nucleate and stabilize micro- 
tubules in their vicinity, a process facilitated by RanGTP, 
which then connect to the astral microtubules emanating 
from the chromosomes, assisted by motor proteins such 
as dyenin (see ,4] for a recent review) . 

Theoretical investigation of the random search-and- 
capture model was done first by Holy and Leibler[5] and 
more recently by WoUman et al[6f]. The latter paper in- 
vestigated the possible effects of a microtubule-stabilizing 
RanGTP gradient around chromosomes, and concluded 
that when many chromosomes are involved, the search 
cannot be completed within the observed time-scales 
without this additional biochemical regulation. However, 
both the papers did not discuss the effects of the cell 
boundary and its constraining effect on the search. Also, 
the analytical results were restricted to the case where 
rescue is completely absent. 

The present paper complements these earlier investi- 
gations into this important problem and extends them in 
several ways. We restrict ourselves to the base search- 
and-capture model (i.e., without considering the addi- 
tional biochemical mechanisms for search), and present 
a complete theoretical solution to the problem of find- 
ing the search time for a fixed target by dynamic micro- 
tubules. Our approach permits the calculation of a com- 
pact mathematical expression in Laplace space for the 
probability distribution function (PDF) of the capture 
time of a single fixed chromosome by successive micro- 
tubules nucleating from a single nucleating site (one at a 
time), for arbitrary cell size and kinetic parameters. In 
particular, we are able to find the search time distribu- 
tion also for non-zero rescue frequency, unlike the earlier 
I^pproaches. We investigate the conditions which opti- 
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mize the search process using this PDF, and conclude 
that the ceh size is an important factor in the optimiza- 
tion process. For small cell sizes, the search time is mini- 
mized when both catastrophe and rescue frequencies are 
small, while for large cells, the rescue frequency need to 
be close to zero, while the catastrophe frequency need to 
be around a certain optimal value. These conclusions are 
found to be in agreement with available experimental re- 
sults for various cell types and existing theoretical results 
where applicable. 

In the following sections, we discuss the problem and 
the model, develop the formalism to address the problem, 
and present our results. The results are then analyzed in 
the context of available experimental observations. 

II. THEORETICAL FORMALISM 
A. The Model 

We approximate the cell as a sphere of radius R, at the 
center of which we place a microtubule nucleating center. 
A single target (chromosome) of cross-sectional area A is 
placed at a distance d from the center. This simplified 
geometry of a mitotic cell permits a complete analytical 
solution of the problem. Microtubules nucleate from the 
center at a rate v. Catastrophe and rescue frequencies 
are denoted by Vc and respectively, and are assumed to 
be the same everywhere inside the cell, as are the growth 
and shrinkage velocities, denoted Vg and Vs- From the 
nucleating center, a single microtubule initially nucleates 
in a random direction and grows by polymerization as 
long as it is in the growth phase, while the same mi- 
crotubule shrinks in length by depolymerization in the 
shrinking phase. In the process, the microtubule scans 
the surrounding space for chromosomes, and is stabilized 
when the growing end encounters a kinetochore. 

A microtubule from a certain nucleation site on the 
centrosome can nucleate in many possible directions; 
however, given the finite size of the centrosome, the ori- 
entation is likely to be constrained by the geometry of 
the centrosome. In an extreme case, one many imagine 
that a microtubule will always grow only along the local 
normal to the surface, but this case is pathological when 
the target is fixed in space, since no microtubule might 
ever grow in the right direction to find it. It is therefore, 
more realistic to imagine that microtubules from each 
nucleation site in the centrosome will grow within a cer- 
tain solid angle Af2, which defines a search cone for the 
corresponding nucleation site. In this case, if the fixed 
target falls inside the cone, and has a cross-sectional area 
a, a microtubule originating at that particular site has a 
probability 

for nucleating in the right direction, within the search 
cone. 



We now need to determine AfJ. We assume that if 
there are M nucleating sites on the surface of the inner 
hemisphere of the centrosome, all the microtubules ema- 
nating these sites together will be searching a solid angle 
47r, without interfering with each other. In this case, we 
require that AfAf2 > 47r, and the inequality allow for 
some overlap of the search-cones in general. To be spe- 
cific, we concentrate on the lower limit of M, and assume 
that M ~ 100 — 300, which fixes the search-cone angle 
as above. In assuming that the search is distributive, our 
model differs from the previous theoretical treatments, 
where each microtubule in the searching population was 
assumed to search the entire solid angle of 47r in the sur- 
rounding space, which, we feel is an unrealistic assump- 
tion. It must also be added that some amount of overlap 
between search-cones may be allowed in our model also, 
but for the time being, we disregard this possibility for 
the sake of simplicity. 




FIG. 1: A schematic illustration of the problem under consid- 
eration. Microtubules nucleate from nucleating sites on the 
centrosome, and search within a certain solid angle about the 
normal to the surface at the site. The single target is as- 
sumed fixed in space. The search cones are assumed to cover 
the entire solid angle, and are assumed non-overlapping. 

A search by a single microtubule is terminated in three 
possible ways: if the nucleation occurs in the right direc- 
tion, (i) the growing tip encounters the target and ends 
the search, or (ii) the microtubule depolymerizes com- 
pletely before it encounters the target, and finally (iii) if 
the nucleation happens in the wrong direction, the mi- 
crotubule will ultimately depolymerize and disappear af- 
ter a futile search, with or without encountering the cell 
boundary. In the last case, we assume that once the mi- 
crotubule hits the cell boundary, it undergoes catastro- 
phe at a rate i'' which is generally higher than the value 
in the interior [3, @, Q . In particular, it was reported in 
[3 that the catastrophe frequency near the boundary is 
16-fold higher than that in the interior, for certain cells. 
In cases (ii) and (iii), a new microtubule will nucleate 
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again from the center in a randomly chosen direction, at 
rate v. We denote by C(T) the probabihty density of the 
capture time, from which the mean capture time follows: 



(T) 



dTC{T)T 
Io°^dTC{T) 



(2) 



In the above equation, dTC{T) gives the probabil- 
ity that the search will be eventually successful. 



B. Capture time distribution 

Since the basic process under consideration here is the 
capture of a certain target by one (or a set) of dynami- 
cally unstable filaments, it is natural to base our th eory 
on consideration of first passage probability densities 
We define three conditional first passage probability den- 
sities (CFPD), which completely describe the search pro- 
cess. These CFPDs are defined below, with their corre- 
sponding conditions in italics. 

1. T) is the CFPD for a freshly nucleated micro- 
tubule to nucleate in the right direction and reach 
a distance d for the first time after a time inter- 
val T, without ever shrinking hack to the origin in 
between. 

2. {\-p)Qx{T) is the CFPD for shrinking to the ori- 
gin after a life-time T, after nucleating in a wrong 



direction, without ever reaching a length X in be- 
tween. 

3. (1 -p)*(r) is the CFPD for return to the origin 
after a time interval T, following nucleation in a 
wrong direction, after encounter with the bound- 
ary (and consequent catastrophe) at least once (and 
possibly several times) in between. 

We now denote by (1 — p)^ftn{T) the probability 
of n unsuccessful nucleation-search-disappearance events 
within a time interval T, prior to a successful event, so 
that 



C(T) = V / n^iT- T')p<i>{d, T')vdT' (3) 



Let us now assume that among the n unsuccessful nu- 
cleation events, there are ni events where the micro- 
tubule nucleated in the right direction, but did not reach 
the chromosome, n2 events where it nucleated in a wrong 
direction, but shrank back to origin before encountering 
the boundary, and n3 = n — rii — n2 events where the 
microtubule nucleated in a wrong direction, encountered 
the boundary, underwent catastrophe and then shrank to 
the origin. 

In terms of the FPDs described above, we find that 
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f]„(r)= / vdTle-"^^ / dT^K^{T^). 
h Jo 



T-Ti'-..T„_i 



(4) 



r 



where K{T) could be Qd{T), QR{t) or '^{T). The tion,a generalized form of the standard convolution the- 
preceding equation has the form of a 2n-fold convolu- orem may be applied here. The result is 

I 



fi„(s) 



1 ( v 



{s + v)\s + v 



{[K,{s)= }^ }^ 



n (■^ + ^) V '5 + ^ 

m —0 n2— ' ^ 



s -\- v I \ s -\- V 



(5) 



From Eql3 we find C{s) = fp^d, s)Y,n=oi^ - 
p)"f2„(s). After using Eq lllBI and performing the ge- 
ometric summation above, we arrive at the following ex- 
pression: 



C[s) 



iyp^{d, s) 



s + iy{l- pQa{s) - (1 - p)[Qr{s) -f 



(6) 
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In particular, if we assume that MT nucleation is suffi- 
ciently fast and not therefore, rate-limiting, we arrive at 
the limiting form 



p^{d, s) 



1 



pQdis) + il~p){ Qfl(s) + vl/(,s) 



(7) 



Eqini is the central result of this paper. Using this 
expression, the mean search time for a single target, and 
its variance may be expressed as 



of appropriate boundary conditions at the origin (nucle- 
ating site) and this has been done earlier (12ij. However, 
we deem this unnecessary for our purpose, since we are 
only interested in using these Green's functions to com- 
pute the CFPDs introduced above. Therefore, for the 
rest of this paper, we will allow the 'length' x to be a 
continuously varying variable between positive and neg- 
ative values, with no boundary condition imposed on the 
dynamics at 2: = 0. The boundary conditions are used in 
the definition of the CFPDs later. 

The Dogterom-Leibler[lli rate equations for MT ki- 
netics takes the form 



(T) 



dsCis)\s=o 
CiO) 



1 d^Cjs) 
(7(0) ds^ 



(8) 



where C(0) is the probability that the search is even- 
tually successful, which we find to be unity in all cases 
studied here. 

The FPDs introduced above are now calculated from 
the Green's functions for MT kinetics, derived explicitly 
in the next section. 



Green's functions 



dtGij — —VgdxGij + VrGoj — VcGij 
dtGoj = VsdxGoj + VcGij - VrGoj (9) 

The equations may be solved together using combined 
Laplace-Fourier transforms, defined as Gij{k, s; xq) 



POO 

J — c 



dx dte ''*Gij{x,t; xo,0). The solution is 



(I,, s) - ^^^^^^2 _ ,kA(s) + Bis)] ^ 



The stochastic state of a MT at a given point in time 
t is characterized by two variables, its length / and its 
state of polymerization versus depolymerization, which 
we denote by an index i, which takes one of the two val- 
ues,! or respectively for growing and shrinking states. 
In this case, therefore, we need to compute four Green's 
functions, or propagators, Gij{x, t; xq^O), for i,j = 0,1; 
by definition, Gij{x, t; xq, 0)dx gives the probability that 
a given MT will have length I between x and x + dx, and 
will be in state i at time t, provided that it had a length 
xq and was at state j at an earlier time t = 0. 

Calculating the above Green's functions for a physi- 
cally realistic situation would also require specification 



where A(s) — [vsi^c—Vgh'r + s{vs—Vg)]/vsVg andB{s) — 

[s{s + Vr+ yc)\/VsVg. 

For connection with the FPDs introduced earlier, it is 
convenient to define the Green's function in such a way 
as that they have dimensions of inverse time, and not 
inverse length. This is done by defining 



= VgGxj ; Fqj = VsGoj (11) 

It will be convenient for later calculations to carry out 
the inversion fc — > a; explicitly: 



Fij{x,s]Xo) 



Foj{x,s;xo) 



Vs [as + (3s 
{as+(3s) 
Vg [as + Ps 
{as + (3s) 



-a.(x-xo)Q(^ - xo) + e^=(^-^°)e(xo - x) 



(12) 



where 

as = ^ + VBis)+A^s)/4 
P.^-^ + VBis)+A^s)/4 (13) 
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and Q{x) is the usual step-function: 8(x) = 1 for a; > 
and otherwise. 



Similarly, 



D. General relations between FPD and Green's 
functions 

The Green's functions calculated in the last section 
may now be used to compute the CFPDs which we used 
before. For this purpose, it is convenient to define first a 
set of unconditional &rst passage densities (denoted FPD) 
as follows: let Cij{x, t; xq,0) denote the probability, per 
unit time, for a MT in state j and with length xq at time 
t = 0, to reach a length x for the first time at time t, and 
in state i. 

For I > 0, Cii (Z, i; 0, 0) is given by the imphcit equation 



Cio(d,s;0) 



Coi{0,s;d) 



-Pad 



(20) 

After solving EqlTT] and Eq[TH] together, we find the 
explicit expressions 



$(d,s) 

Qd{s) = ■ 



D{s)e- 



VrVc[\ - e-("=+/3=)'i] + D{s) 



1 - e 



(21) 



where -D(s) is given by 



Fn{l,t;0,0) = Cn{l,t;0,0) + 
lim / dt'Cii{l,t';0,0)Fii{l~€,t;l,t') (14) 

In the above equation (and the following equations), 
the e-factors take into account the following restriction 
on its dynamics: starting from a growing state at i = 0, 
with a length I, it can return to the same length I at a 
later time, in a growing state, only from below (which 
decides which of the 8-functions appearing in Eq[T2] is 
non-zero). Similar restrictions apply to the equations 
below. 

Similarly, Coi(0,T;d,0) and Cio(d, T; 0, 0) are given 
by the equations 



D{s) = {s + asVs)is + PsVg) + l^r{s + (3sVg) + 1^c{s + asVs)- 

(22) 



E. Catastrophes at the cell boundary 

We assume that when a MT hits the cell boundary 
by growing, it undergoes catastrophe there at a rate v'^. 
We now compute ^'(r), which is the CFPD of return to 
origin (i.e., complete depolymerization) of a MT after a 
lifetime T, and an encounter with the boundary at least 
once. 

Clearly, along the line of our previous arguments, ^'(r) 
may be given by the expression 



Foi(0,T;d,0) = Coi(0,T;d,0) + 

lim / dT'Coi{0,T';d,0)Foo{e,T;0,T') (15) 

Fio(d,T;0,0)-Cio(d,T;0,0)-H 

lim / dT'Cio{d,T';0,0)Fn{d-€,T;d,T') (16) 

Using these two FPDs, we are now in a position to 
write down the following relations between the CFPDs 
introduced earlier: 



*(T) = [ dTi^{R,Tx) X 
Jo 

fT-Ti 

/ v'^dT2e-'''^^-x{R,T-Ti-T2), 
Jo 



(23) 



where x(^? T) gives the FPD of complete depolymer- 
ization of a MT, starting at the boundary, at length R 
in shrinking state, with possibly multiple visits back to 
the boundary in between. This quantity may now be 
expressed implicitly through the equation 



Cii(d,t;0,0) = $(d,0+ / dt'Qd{Q,t')Cio{d,t-Q,t') 

Jo 

(17) 



Coi(0,T;0,0) = Qd(r)-f / dr'$(d, r')Coi(0, T; d, T') 

Jo 

(18) 

Eg [MlfTSl mav now be solved using Laplace transforms. 
From EqlHl we find that 



Cii{d, s; 0) = lim 



Fii(d,s;0) 



£-0+ 1 + Fii{d~ e,s;d) 



(19) 



X{R,T) = ^*{R,T) + [ dTiQ*j,{Ti) X 
Jo 



T-Ti 



i^',dT2e-'''^^^x{R, T-Ti~ T2) (24) 



where Q*j^{T) is a 'mirror' image of Qr(T) introduced 
earlier, and represents the CPFD of a return to boundary 
over a time interval T, without ever reaching the origin 
(i.e., shrinking to zero) in between. Similarly, $*(i?, T) 
is the 'inverse' of <f>(i?, T), and gives the CFPD of com- 
plete depolymerization of a MT starting at the boundary, 
without ever returning to the boundary in between. 
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Eq|23l and Eq[24] may now be solved together using 
Laplace transforms, and we find 



Quantity Yeast (I) 



Mammalian(II) Xenopus{lll) 



iy'^^{R,s)^*{R,s) 



1 - QU-') 



(25) 



l>*(i?, s) and (9^(s) may be obtained from s) and 



Qr{s) respectively by the transformations Vs 



and 



— > (and vice-versa). From Eq[T3l this has the effect 
of replacing as by /3s and vice-versa, while D(s) defined 
in Eq l22l remains invariant under these transformations. 
Therefore, using Eql2H we arrive at the following expres- 
sions for these 'inverse' quantities: 



asVs 



$*(i?, s) 



1 - e 



D{s)e 



-PsR 



UrVc[l - e-("=+/5.)«] + D{s) 



III. RESULTS 



(26) 
(27) 



The mean time of capture of the target is now evalu- 
ated using Eql7|and Eq [8] directly in Laplace space, since 
explicit inversion of Eq[7] turned out to be too difficult, 
and an explicit expression for the mean time in the most 
general case turned out to be too lengthy for reproduc- 
tion here. We, therefore, used MATHEMATICA (Wol- 
fram Research Inc., version 7) to compute the mean time 
from EqlH after substituting Ea l2Tll22l and Eq |25ll27l In 
all the following analysis, we use the limits v ^ oo and 
u'^ — !■ cx), i.e., microtubule nucleation is assumed very fast 
and boundary-induced catastrophes are assumed to be 
instantaneous. These assumptions are not strictly neces- 
sary, but reduces the number of active parameters. Since 
a finite nucleation rate and a non-zero delay before catas- 
trophe at the boundary would increase the total time, our 
estimates for the mean search time obtained henceforth 
may be regarded as theoretical lower limits. 

Experimental measurements of the microtubule kinetic 
parameters show some distinct c hang es as the cell pro- 
gresses from interphase to mitosis [Ijl. However, the na- 
ture of the changes are different for different types of 
cells, and the general qualitative features appear to bear 
a strong correlation to the cell size. Yeast cells show 
a reduction in both catastrophe and rescue frequencies, 
but the changes are relatively small. In mammalian cells, 
which are typically larger, there is a marked fall in rescue 
frequency between interphase and mitosis, and a two-fold 
increase in catastrophe frequency. In Xenopus oocytes 
(frog egg cells, which are large and almost a mm in ra- 
dius), the effects are somewhat different; while the res- 
cue frequency, while small, is almost doubled, there is a 
sharp, seven-fold rise in the catastrophe frequency. 

In order to understand the optimization criteria used 
by cells, we investigated the mean search time for three 



Ur{exp) (0.42)0.12 min-^ (10.5)2.7 min ^ (0.66)1.2 min ^ 

(exp) (0.48)0.24 min"^ ((1.56)3.48 min"^ (1.08)7.2 min-^ 

Vg (exp) 1.7 j-im min^^ 12.8 /im min^^ 12.3 fim min~^ 

Vs (exp) 2.7 /im min~^ 14.1 /im min~^ 15.3 /im min~^ 

R 5/im 20/im 100/im 

d 10/im 20/im 

TABLE I: Experimental values of microtubule kinetics in the 
mitotic phase, as taken from [l3l].The values in brackets are 
those in the interphase, prior to entering mitosis. The values 
for the cell size and distance to the target (chromosome) are 
only rough estimates. For theoretical analysis, we used Vg = 
2.0 and Vs — 3.0 /im min~^ for I, and Vg = 12.0 and Vs = 14.0 
/im min~^ for II and III. 



sets of parameter values, corresponding to Yeast (I), 
mammalian(II) and Xenopus(III) cells, whose experi- 
mentally observed values are summarized in Table I. Of 
these, we use only the growth and shortening velocities 
precisely, while Vr and Vc are varied in the range 0-5/im 
min~^. The cell sizes given are only approximate es- 
timates. For case III, we used a somewhat smaller cell 
radius to avoid over-flow in computation. In addition, we 
also study a theoretical case where i? — > oo separately, 
for which an exact expression for the mean search time 
is obtained in Appendix I, when — 0. 



A. Contour plots of mean search time 

FiglH displays a contour plot of the mean search time 
for case I, and Figl3] gives the same for case II, with M = 
100 searching microtubules in both cases. We verified 
that when the number of microtubules is doubled, the 
search time is almost halved in all these cases. It may 
also be noted that our contour plots are similar with those 
obtained numerically in 



B. Optimization of Uc at fixed i/r 

We now check how the search time varies as a function 
of Vc, for various fixed values of Vr, i-e., horizontal cross- 
sections of the contour plots presented earlier. 



1. Small cells: R — 5/im, d — 3^j,m 

FigU shows the variation of the mean search time ver- 
sus catastrophe frequency for three different Vr, corre- 
sponding to Figl2] earlier. The smallest search times are 
obtained for Vr — and Vc — 0, implying that for small 
cells, the search is optimized by ensuring that the micro- 
tubules simply grow all the way to the cell boundary and 
undergo catastrophe, but with no rescue or catastrophe 
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mean time(inin) 
35 ■ 





v,=0.0 




v,=0.25 




v,=0.5 




v,(l/min) 



yt-(l/min) 



FIG. 4: The mean time of capture of a single chromosome by 
M = 100 microtubules for a cell of radius R = 5/im. The 
chromosome is fixed at a radial distance of d = 3^tm from the 
centrosome. 



FIG. 2: Contour plot of the mean time of capture (in mirmtes) 
of a single chromosome hj M = 100 microtubules for a cell 
of radius R — 5fim. The chromosome is fixed at a radial 
distance of d = 3/im from the centrosome. 




12 3 4 

yj(l/iiiin) 



the mean search time is observed to occur at the extremal 
values Vr — O^Vc = 0, while for Vr > 0, the minimum 
occurs at a non-zero value of Vc, which increases with v^. 
Fig|S] shows similar curves for a larger cell: R = 20/im 
and d = 10/im, and with more searching microtubules; 
M — 400 in this case. The number of microtubules were 
chosen such that the minimum of the mean search time 
remains less than 30 minutes. 

When comparing these predictions with the experi- 
mental data in Table I, we observe that, while i^c doubles 
between interphase and mitosis, Vr is reduced by almost 
a factor of 4. This may be interpreted as follows: since Vr 
is high in the interphase, it is more important to reduce 
it to levels low enough that the search does not become 
prohibitively costly, time-wise, and the relatively smaller 
change in Vc is likely to be a fine-tuning exercise. 



FIG. 3: Contour plot of the mean time of capture (in minutes) 
of a single chromosome by M = 100 microtubules for a cell 
of radius R = lOnm. The chromosome is fixed at a radial 
distance of d = 5/im from the centrosome. 



in the cytoplasm itself. These observations are in agree- 
ment with the experimental observations in Yeast cells 
(Table I), where both these events are seen to become 
rarer in mitosis as compared to interphase. 



mean timc(inin) 
30 





v,=().0 




v,=1.0 




v,=2.0 




[■..(l/min) 



2. Large cells:R = 10 — 20/im, d = 5 — lOfim 

We next consider somewhat larger cells of radius 10 — 
20/im. Figl5]shows a cross-section of Figl3]at three rescue 
frequencies. In this case, again, the absolute minimum in 



FIG. 5: The mean time of capture of a single chromosome by 
M = 100 microtubules for a cell of radius R — lOpim. The 
chromosome is fixed at a radial distance of d = 5/im from the 
centrosome. 
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FIG. 6: The mean time of capture of a single chromosome by 
M = 400 microtubules for a cell of radius R — 20^m. The 
chromosome is fixed at a radial distance of d = 10pm from 
the centrosome. 



3. Very large cells:R — lOO^m, d — 20fj.m 

As the last case, we consider very large cells, with ra- 
dius R =^ 100/im. This case is motivated by Xenopus egg 
cells, which are known to have radii of the order of 1mm. 
We fixed the chromosome separation d — 50/im here, 
and increased the number of searching microtubules to 
M = 1000. FiglT] shows the results. In this case, except 
for Vr — 0, we see a pronounced minimum of the mean 
search time as a function of Vc- This trend agrees with 
the observations in Xenopus extracts (Table I), which 
show a modest increase in rescue frequency, but an al- 
most seven- fold increase in Vc- We interpret the latter 
as arising from the requirement to tune i^c in order to 
minimize the search time. Note also from the figure that 
the curve is very strongly asymmetric on either side of 
the minimum: a larger Vc would increase the search time, 
but to a much lesser extent compared to a lower Vc- 

We also studied the fluctuations in the ti me of capture , 
by looking at the RMS deviation AT = \/{T'^) - (T}2. 
FigEK shows the relative fluctuation, i.e., AT/(T) for 
the three different cases of small, large and very large 
cells studied here, when i^r = 0. The relative fluctua- 
tion stays close to 1 for the entire range of studied, 
but when the rescue is increased to i^r = l.Omin"^, the 
relative fluctuation increases by 10 to 20% (Figl5}3. The 
trends in the relative fluctuation indicates that, at least 
for small i/^, the probability distribution of the search 
time might be exponentially decaying (a distribution of 
the form C{T) = Xe^^'^ has a mean of and relative 
fluctuation 1); however, more investigations are needed 
here. 



IV. CONCLUSIONS 

In this paper, we introduced a rigorous first-passage- 
time-based stochastic formalism to study the problem of 
calculating the time of capture of a fixed target by a set of 



mean time(min) 
700 [■ \ 



^,.(l/mln) 



B. 



FIG. 7: (A). The mean time of capture of a single chromosome 
by M = 10^ microtubules for a cell of radius R — lOO/im. The 
chromosome is fixed at a radial distance of d = 20pm from 
the centrosome. (B)The mean time for Vr = 1.0, for a smaller 
range of Vc- 



dynamic microtubules, nucleating from a center (the cen- 
trosome). In contrast with earlier theoretical approaches, 
we have been able to obtain estimates of the mean search 
time for arbitrary rescue frequency and cell size. 

Our results show that the cell size is an important fac- 
tor in determining the optimization of microtubule kinet- 
ics in the mitotic phase so as to minimize the search time. 
For very small cells, with size R < 5/zm, it is advanta- 
geous to minimize both the cytoplasmic catastrophe and 
rescue events, with only the boundary terminating the 
search. In larger cells, the primary condition is to mini- 
mize cytoplasmic rescue, and a certain amount of tuning 
of Vc is advantageous, though it is not catastrophic if the 
actual Vc is somewhat smaller or larger, as the minimum 
is quite shallow. In much larger cells, the search time 
becomes a rather sensitive function of Vc, and tuning Vc 
becomes the priority; while Vr is required to be as low 
as possible, changes in Vr does not affect the search time 
as much as changes in Vc- In particular, the search time 
rises rather steeply below the optimum value. The trends 
in microtubule kinetics as predicted by this analysis are 
in good qualitative agreement with the existing experi- 
mental results. 

We are presently working on extracting information 
about the probability distribution C (T) from our formal- 
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FIG. 8: (A). The relative fluctuation in the time of capture 
of the target for three cases, I, II and III (Table I), when 
Vr = O.Omin"^(A) and Vr = 1.0min"^(B). 



ism, and extending our analysis for multiple targets, and 
hope to report our results in the near future. Other pos- 
sible extensions of this study involve including target dif- 
fusion and side-capture of microtubules by chromosomes 
via motor proteins (distinct from the end-capture mecha- 
nism explored in this and previous theoretical models). 
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APPENDIX A: CALCULATION OF MEAN 
SEARCH TIME WHEN Vr = Q AND R = oo 

When i/r = 0, we find that 

[s + Vc)/vg ; P^s/vs, (Al) 
from Eq[T31 Upon substitution into Eql221 this gives 



D{s) = + e's + 6, 



(A2) 



where 



e = vl^ ■ e' ^2v,{l + vslvg). (A3) 

Let us now consider Eql?] When R ^ oo, ^{s) — 0, 
and Qr{s) — > Qoo(s), where 



Qoo(s) 



Vc + s(l + Vg/Vs) 



(A4) 



from EqEIland using Eq|A2land EqEHl above. There- 
fore, we have the following limiting form for C(s): 
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Coois) = p- 



1 - 



pQd{s) + (1 -p)Qoo(s) 



(A5) 



If we now expand Coo (s) = Cq — sCi + ..., it is clear that 
(T) ~ Ci/Co from EqEl We can therefore substitute 
Eq lATlAil in EqEH substitute into EqEll and expand 
the resulting expression upto the linear term in s. The 
result is Co = 1 and Ci = —{T)ao, where 



in [01 , using a somewhat less rigorous method. In this 
limit, the search is optimized at i>c = d/vg. 

The expression in Eq lA6l above may be approximately 
reproduced by physical arguments as follows. The prob- 
ability that a microtubule will nucleate in the right di- 
rection, and will not undergo catastrophe until it reaches 
the target is given by Ps = pe~'^'''^^^<' , and it will take 
at least N ^ pj^ unsuccessful attempts before this is 
accomplished. Each of these unsuccessful search events 
lasts a time r , and therefore, the total search time 

is 



(T)c 







re". 1 










Vs) 


p 



d 



(A6) 



is the mean search time in an infinitely large cell, with 
no rescue. Note that, in the limit p — > 0, the above 
expression reduces to (1 -I- Vg/vs)e^'='^^^<' /pvc, as derived 



PVc 



(A7) 



Note that Eq lA6l reduces to Eq lA7l in the limit p — > 



and Vg < Va 



